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We introduce a model of interacting Random Walk, whose hopping amplitude depends on the 
number of walkers/particles on the link. The mesoscopic counterpart of such a microscopic dynamics 
is a diffusing system whose diffusivity depends on the particle density. A non-equilibrium stationary 
flux can be induced by suitable boundary conditions, and we show indeed that it is mesoscopically 
described by a Fourier equation with a density dependent diffusivity. A simple mean-field description 
predicts a critical diffusivity if the hopping amplitude vanishes for a certain walker density. Actually, 
we evidence that, even if the density equals this pseudo-critical value, the system does not present 
any criticality but only a dynamical slowing down. This property is confirmed by the fact that, 
in spite of interaction, the particle distribution at equilibrium is simply described in terms of a 
product of Poissonians. For mesoscopic systems with a stationary flux, a very effect of interaction 
among particles consists in the amplification of fluctuations, which is especially relevant close to the 
pseudo-critical density. This agrees with analogous results obtained for Ising models, clarifying that 
larger fluctuations are induced by the dynamical slowing down and not by a genuine criticality. The 
consistency of this amplification effect with altered coloured noise in time series is also proved. 

PACS numbers: 05.60-k Transport processes - 05.40.Fb Random walks and Levy flights - 05. 40. a Fluctuation 
phenomena, random processes, noise, and Brownian motion 



I. INTRODUCTION 

Diffusive dynamics is at the core of the study of heat 
or matter transport 0, In ia,ct, exceptions to this 
paradigm are themselves matter of research Q. Typ- 
ically, one expects that non-equilibrium steady states 
are described by a Fourier-like equation with non con- 
stant diffusivity; for homogeneous systems, such a dif- 
fusivity should depend on the local density of energy or 
matter. Actually, this kind of behaviour was revealed 
in different fields - magnetic systems, biology, reaction- 
diffusion processes (see for example [1, H, [a, Q)- On a 
more modellistic ground, we may quote a recently stud- 
ied 2-dimensional Ising system between two thermostats, 
where heat transport exhibits indeed an energy depen- 
dent diffusivity [i, i, 0, nil • 

A basic question is the relation between the mesoscopic 
description given by the Fourier equation and the under- 
lying microscopic dynamics. In this context, geometry, 
interactions and possible criticalities should also be taken 
into account. In fact, various features depend on this 
mesoscopic-microscopic relation, e.g. scaling properties, 
behaviour of fluctuations, coexistence of critical and non- 
critical phases within steady states, etc. 

In order to carry out the analysis, a suitable resource 
is given by the prototype of diffusive discrete processes, 
i.e. random walks (RW) [l^l- In its original form, a par- 
ticle (the walker) has a constant uniform probability to 
jump at discrete times from a site to another in its neigh- 
borhood. The addition of constraints (e.g. self avoiding 
rules) or the dependence of the probability on the site or 



the neighborhood state can fit a variety of interesting sit- 
uations, including non-equilibrium dynamics. Moreover, 
particle reservoirs with different densities may be fixed 
at borders, determining a fiux from the "hot" (dense) 
to the "cold" (rare) side. Therefore, we expect that in- 
teracting RW could represent various kinds of transport 
properties. In this way, the particle population can be 
read as a diffusing "energy" . Such an interpretation fits 
the behaviour of many diffusive systems, physically very 
different in principle, from sandpiles to spin models, as 
those quoted above. 

A number of conceptual achievements have been ob- 
tained in these fields focusing in particular on the sym- 
metric simple exclusion process (SSEP), where two walk- 
ers cannot occupy the same site, and on the zero-range 
processes (ZRP) where the hop rate depends on the parti- 
cle density of the departure site (see reviews [H, Q)- An 
obvious extension of ZRP are the misanthrope processes 
(MP), where the hop rate depends also on the density 
of the target site. Actually, MP have been investigated 
both for the factorization conditions of equilibrium states 
and for the mapping onto other models. Moreover, for 
SSEP and ZRP the situation where a stationary flow is 
imposed by the presence of particle reservoirs has also 
been considered. In particular, the problems of the sta- 
tionary density profile 15] and of current fluctuations 
[l6; 17, 18] have been deeply investigated. 

Here, we want to focus on the possibility of modelling 
some of the features characterizing a diffusive dynamics 
with non constant diffusivity. We would give evidence, 
for instance, to the relevance of a possible critical phase in 
mesoscopic phenomena, especially those regarding fluctu- 
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ations and decay of correlations. 

Consider for definiteness a one dimensional discrete ar- 
ray with particle reservoirs at the extremes (the general- 
ization to higher dimensional lattices is straightforward) . 
We shall introduce a particular case of symmetric MP 
where the hop probability is "totalistic" : this means that 
it depends only on the mean value, hereafter called "link- 
population" , between the departure-arrival populations. 
Such a model is expected to describe a general diffusive 
system where the local conductibility of a microscopic 
connection depends on the average density of the link it- 
self, i.e where diffusivity is density dependent. Moreover, 
for a possible critical behaviour, the dynamics should also 
include a parameter Nc, such that the hop rate tends 
to as the link population approaches Nc- In a mean 
field description, this implies a vanishing diffusivity at 
Nc, suggesting indeed a possible criticality. 

Analytical and numerical investigations show that our 
RW system actually exhibits diffusion, being well de- 
scribed by a Fourier equation for the population density. 
In particular, in agreement with our intent, the diffusiv- 
ity depends on the average population. Moreover, there 
is a dynamical slowing down when this population equals 
Nc, and mesoscopic realizations prove that there exists a 
band around this putative critical parameter where fluc- 
tuations are amplified with respect to the equilibrium 
values. However, it is possible to show that Nc is not a 
true critical point; for large sizes, in particular, fluctua- 
tions tend to vanish in a clear analytical way. 

This constitutes an interesting result, especially in 
comparison with the spin system mentioned above [Tol . 
[Tlj . There are relevant analogies: the spin flip in Ising 
and the jump probability in RW both depend on " link 
properties " (energy and link populations, respectively); 
both systems display a pseudo-critical parameter (the 
critical temperature or energy, say Ec, and the popula- 
tion Nc, respectively). Nonetheless, since evolution rules 
are different, there are also important differences: not 
only Ec and Nc have a different critical nature, but also 
dimensionality, for instance, plays a quite different role. 
Since for both systems we evidence amplified fluctuations 
at the mesoscopic scale, we can conclude that their source 
is not the (virtual) presence of a genuine critical point, 
but rather the slowing down in dynamics induced by par- 
ticle interaction. 

After the introduction of model and formalism (section 
HI)) , we investigate the system first by analytical methods 
and then by numerical simulations. The latter corrob- 
orate and expand analytical results. More precisely, in 
section Hill we prove that, for a closed system, a hopping 
probability depending on the link population is sufficient 
to derive a factorized Poissonian distribution describing 
the statistics of the local population. Then, in Section 
HVl we take advantage of this result to get an analytic ex- 
pression for the diffusivity, which is valid for sufficiently 
large systems. Such an expression is confirmed by nu- 
merical simulations displayed in Section |Vl We also dis- 
cuss data relevant to fluctuations for proper observables 



and their possible correlations with the quasi-critical re- 
gion, the meaning of such a dynamical slowing down and 
the mesoscopic character of the phenomenon. Finally, 
we study internal correlations by means of the spectral 
analysis of time series. 



II. THE MODEL 

A random walk is a stochastic process defined on a 
generic graph. At each discrete time-step, the walker 
moves to an adjacent site following a random rule which 
is independent of the walk history, so that the process 
is Markovian. For the simple random walk, the walker 
jumps on one of its neighbouring sites with uniform prob- 
ability. 

In a system of N simple walkers (particles) on an Eu- 
clidean lattice the limit of continuous time and space is 
described by the Fourier-like diffusion equation 

dttiix,t) = -l) ■{Dl)iiix,t)), (1) 

where /.t(x, t) is the particle density in x at time t and 
the constant D is the diffusivity. Note that the quantity 
which is derived on the right side of i.e. 

<j>ix)^-iD'^^,{^,t)), (2) 

is the local flux at x. 

Equation ([T]) has been introduced in different contexts 
to study matter or energy transport. In particular, even 
for several interacting models, transport is expected to 
be described by an analogous Fourier-like equation with 
non constant diffusivity [1, 2]. In this work, we aim to in- 
troduce a simple microscopic model where, at mesoscopic 
level, the diffusivity depends on the particle density fi in 
a tunable way. 

More precisely, we consider a system of walkers on a 
one dimensional lattice of L sites i = l,...,L (anyway, 
most of our results can be generalized straightforwardly 
to higher dimensional lattices). We study both the close 
chain with periodic boundaries and the system with open 
boundaries in contact with particle reservoirs. The tran- 
sition probability between neighbouring sites is a sym- 
metric function of and rii+i, i.e. the number of walkers 
occupying sites i and i -|- 1 respectively. In the following, 
Ui will be called occupation number or site energy, and 
Ei = {rii + rii_|_i)/2 will be called link population or link 
energy. In order to simplify formulas, we shall often use 
fii instead of the more explicit (rii). 

Now, let f{ni,ni+i) be the probability for a particle 
on the site i to jump on the adjacent site i + 1. 

In a mean field approach, we discard the number fluc- 
tuations, assuming that Ui ~ {rii). The average number 
of particles crossing the link i,i + I in the unit time, i.e. 
the particle flux, is 

<^{i,i + l) = {n^)f{{n,), (nj+i)) - {n^+i) J {{n^+i) , {n,)). 
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If the transition probability f{ni,ni+i) is symmetric in 
its arguments, we can write: 



d{n) 



dx 



i+l) 



(3) 



where d{n)/dx, or equivalently dfi/dx, represents the 
discrete spatial derivative for the average number of par- 
ticles. 

The previous expression simply states that the flow is 
proportional to the particle gradient. Therefore, it is a 
discrete one dimensional version of the flux ([2]) appearing 
in the Fourier equation ([1]). Now, if /(ni+i,ni) depends 
only on the link energy Ei, and if the system is large 
enough so that the energy differences between neighbour- 
ing sites can be neglected, by comparing Eq.[3]and Eq.[Tl 
we can express the diffusivity as D{fi) — f{^J,). 

It must be underlined that the approach followed in 
this section is a mean-field type, since fluctuations in 
particle density have not been taken into account. In 
the next sections we make a step forward by introducing 
a probability distribution for the particle configurations, 
and we obtain a more reliable expression for the diffusiv- 
ity. 

Finally, we notice that if Eq. [T] holds, the diffusivity 
D{fj,) can be calculated directly from the spatial distri- 
bution of the average density /ii. In fact, in a stationary 
system, the average flux is constant along the direction 
of the underlying structure. Hence we obtain 



(4) 



where the constant C is the flux. 

If for different system parameters (size, boundary con- 
ditions,...) Equation (j4]) provides the same expression of 
the diffusivity as a function of the density, this would be 
a reliable check of the validity of the Fourier equation ([T|) 
in describing matter transport in the system [10|. 



III. ANALYTICAL RESULTS ON CLOSE 
SYSTEMS 

Consider a close system of N particles with the follow- 
ing updating rule: at each time step, a randomly chosen 
particle can move from its site i, towards sites i -I- 1 or 
I — 1 , or else stop, according to the transition probabili- 
ties fr = f{ni,ni+i), fi = f{ni,ni-.i) and fs = l-fr-fi, 
respectively. Hence, when the particle moves from i to, 
say, i + l, then rii changes into = — 1 while rii^i 
changes into n'^^^i — n^+i + 1. 

The stationary state of such a system is described by 
the distribution P({n}), where {n} represents the parti- 
cle configuration defined by the occupation numbers of 
each site. Obviously, ^ ■ rii = N. 

The detailed balance condition reads: 



(5) 



where ^ represents the probability that one of the parti- 
cles on i is selected. The solution of this equation would 
provide the probability distribution for the equilibrium 
system. 

As already noticed, if a transition probability is to- 
talistic, i.e. it depends on the sum rii + then 
f{n„ni+i) = In fact, n^+j + nj; = Ui+i + 

1 + ni — 1. Therefore, Eq. Ogets 



P{{n})n, = 



(6) 



We note that totalistic transition probabilities satisfy 
the condition given in reference [isj for the factorizability 
of the distribution function in a MP. Therefore we can 
assume: 



P{{n})^Y[P,{r 



(7) 



The previous equation allows to focus, in Eq. [SI on the 
factors relevant to indices i and i + l only, giving 

By properly separating terms we obtain the following re- 
cursive equation: 



P^{ni) 



P.{n, ~ I) 



Pi+i{n 



i+l 



1) 



Pi+i{ni+i) 



Tii+i — const, (9) 



whose solution is just a Poissonian distribution indepen- 
dent of the site index (the constant does not depend nei- 
ther on the site nor on the particle number). More pre- 
cisely, we have 



P^{n^)=V^,{ni) 



(10) 



where fi is the average number of particles per site. 

Clearly, in the presence of a density gradient, the pre- 
vious expression is no longer correct. In fact, the mean 
occupation number /i is in principle site dependent. How- 
ever, we expect that the Poissonian distribution still pro- 
vides a good approximation when the system is suffi- 
ciently large. 

Finally, it is important to notice that the request im- 
posed on the transition probability / (i.e. to be depen- 
dent on the sum -t-n^+i only) proved to be sufficient to 
get the distribution pO)l . We remark also that the equi- 
librium distribution depends on the particle density but 
not on the form of /(n^+i, rii). The Poissonian distribu- 
tion therefore characterizes a wide class of random walk 
systems, and not only the simple non-interacting model. 



IV. ANALYTICAL APPROACH TO 
DIFFUSIVITY 

Let us now consider a system endowed with a density 
gradient: two particle reservoirs are placed at the bound- 
aries forcing the occupation numbers of border sites to 
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assume different values. More precisely, at each time step 
the particle numbers on sites i = 1 and i = L are drawn 
from Poissonian distributions with averages Ni and Nl 
respectively (iVi < Nl). 

The flux between sites i and i + 1 (Eq. [3|) has to be 
rewritten taking into account the fluctuations of the par- 
ticle numbers: 

+ = ^P{{n})n^f{n, + ni+i) - 
{»} 

-^P(MK+i/(n,+n,+i) (11) 

{n} 

The previous expression is exact, since it takes into ac- 
count all possible configurations. Conversely, Eq. [3l 
where occupation numbers are treated as average values, 
is a mean-field approximation of Eq. 1111 

Now, in the presence of a density gradient, the analyt- 
ical expression for the configuration probability P{{n}) 
is not known. However, if we consider a large system 
such that Nl ~ Ni ^ L, a, reasonable ansatz is to as- 
sume that the distribution is still factorizable and that 
in each site the occupation probability is described by a 
site dependent Poissonian. 

The average occupation number of site i will be provi- 
sionally denoted fi instead of fii. Assuming a certain 
smoothness, for the adjacent site we can consistently 
write fJ.i+1 — IJ. + s, with e small. 

Within such approximations, the first term in (jlip is 

nf{n + m)'Pf,+e{m) = 

e-(A'+'^)(^ + e)™ 



E 

n,m 

E 



n! 



71/(71 -|- m)- 



m! 



71/ (n -I- m) 



(1 - e) 1 + -771 







1 + 









(12) 



where we neglected orders e^. 

The second term in Eq. [Tl] can be rewritten analo- 
gously as: 



^'Pfj.+eim) ■mf{n + m)Vf,{n) 
y^V^{m)V^{n)mf{n + m) 

n.ra 

By subtracting Eq. [H and Eq. [T3l 

^7',,(77i)7'^(7i)/(71 + to)- ( 







1 + 









77171 — 771 



13) 



(14) 



Notice that, due the symmetric sum, we can write 



mn — m as — -^^ — t^-^ 



. Therefore we obtain: 



eE p^/(n + 7n)— (m-7i)' = $(i,i + l), 

^ — ' rj!77i! '■■ 



2fi 



(15) 

to be compared with Eq. [31 The analogy is actually 
straight as e is just dn/dx. Interestingly, from Eq. [15] 
we obtain that the approximate diffusivity D is 

~ , „ e~2p n+m rm-7l)^ 



2/i 



confirming that D only depends on ^. However, we recall 
that this is still a local relation, since we posed = /i. 

Provided that ?7i and 71 are sufficiently large, a sim- 
ple analytic relationship between D{^) and f{n -t- tti) is 
obtained approximating the Poissonian by a Gaussian 

distribution {V^{n) ~ ^ ) and adopting a 

continuous picture: 



00 /'OO 



(17) 

/(a; y)dxdy. 



/2¥]2 y'27r/i 2/1 
Now, the variable change: a; 4- y = and x — y — t] gives: 



^-00 

/"OO 



87r/Lt2 

(18) 



It is now suitable to explicit the transition probability 
/(Tii, 7ii+i). We recall that it should depend on the link 
energy and vanish for Ei = Nc . We therefore define 



1 r 
2 



1 



(19) 



and this expression will be used in numerical simulations. 
Once inserted in Eq. I18[ such a transition probability 
provides a close form for the diffusivity: 



1 



(20) 



Notice that D{fi) is qualitatively analogous to /(/i) (Eq. 
[T9| but it does not vanish for fi = Nc- Consequently, 
1^ does not diverge and Nc is not a true critical density, 
but only a value where the dynamics presents a slowing 
down. The absence of critical phenomena was already 
foreseeable by considering that at the equilibrium the 
distribution is regular for any value of fj, and Nc- 
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V. NUMERICAL RESULTS 



Before discussing the results obtained from simula- 
tions, we briefly resume how numerical experiments work. 
A system, characterized by parameters {L, Ni, N^, Nc), 
is prepared with a uniform distribution of particles. At 
each unit time, a randomly extracted walker hops accord- 
ing to the transition probability given by Eq. 1191 Then, 
the system evolves along this rule, generating time series 
for each observable. Such series evidence how the system 
looses the memory of the initial distribution, eventually 
reaching a steady state that depends only on the system 
parameters. To any observable we can associate its mean 
value (•) and fluctuation F(-) . Such quantities are eval- 
uated by considering a large number (~ 10^) of different 
independent realizations. 

A. Diffusivity 

Let us first verify that the system can be described by 
means of a Fourier like equation with a density dependent 
diffusivity. For this purpose we show that equation Q 
holds for different boundary conditions and system sizes. 
In Fig. [T] we have plotted the value of the discrete deriva- 
tive {d{n) /dx)~^ (multiplied by a suitable constant) as 
a function of (n) . Different markers correspond to dif- 
ferent boundary conditions and sizes. The data collapse 
is very good, confirming the general picture of a density 
dependent diffusivity. The continuous line represents the 
theoretical prediction given by Eq. [501 We recall that 
Eq . EOl was derived under the assumption that the global 
distribution P{{n}) is factorized into L Poissonians Vfj.^, 
moreover we approximated such Poissonians with Gaus- 
sian distributions. Within such approximations, the non 
perfect coincidence between the analytic and the exper- 
imental minima can be considered as natural. The cal- 
culated expression D{n) provides therefore a good esti- 
mate of the diffusivity, even for small sizes and occupa- 
tion numbers. In Figure [2l such a good agreement is 
evidenced for different values of Nc. 

Looking at the well collapsed minima, numerical ex- 
periments confirm that the diffusivity is strictly positive: 
consequently, we do not expect any critical effects. More 
precisely, the site ic (i.e. the site approximately corre- 
sponding to the density Nc) does not give rise to any sin- 
gularities in the observables describing the system. Fur- 
ther evidences are given in the next sections. 

B. Poissonian 

We study now the probability distribution for the par- 
ticle population in the presence of a gradient. As proved 
in Section lTTTl for a closed system the occupation number 
is a stochastic variable described by a Poissonian distri- 
bution. For an open system with a steady particle flow, 
the average occupation number is site dependent and its 




FIG. 1: Simulations for Nc = 8. Data points (symbols) and 
theoretical results D{{n)) (line) for C{d(n) / dx)~^ \ different 
symbols refer to systems of different sizes and boundary con- 
ditions, as explained by the legend. The constant C has been 
fixed for different data in order to obtain the curve collapse. 
The theoretical curve is given by Eq. [20] 
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FIG. 2: Data points (symbols) and theoretical predictions 
(lines) for the diffusivity as a function of the particle density 
for two values of the parameter N^- The numerical data are 
obtained for systems of size L = 64, with A^^i = 2 and = 10. 



distribution is no more rigorously described by such a 
Poissonian. However, we expect that the differences be- 
tween the exact distribution P({n}) and small 
when the size L is large with respect to the particle gra- 
dient Nl — N\ . We also expect that, the farthest from Zc, 
the better the approximation works. Indeed, these facts 
are confirmed by numerical simulations. For example, in 
Fig. [3| we show the probability distributions for the oc- 
cupation number relevant to several sites of an open sys- 
tem of size L = 16, with iVi = 2, TVl = 10 and N^ = 8. 
Each site i considered (depicted in a different colour) cor- 
responds to an average occupation number /i^, and the 
pertaining distribution is compared with the Poissonian 
of average /i^, i.e. the distribution for a close system made 
up of iV = L X /ii particles. As expected, near the bor- 
ders the two distributions perfectly overlap; in general, 
their difference is small away from (which is about 12, 



FIG. 3: Occupation number distributions P{i, k) for L = 16, 
A^^i = 2, Nl = 10 and Nc = 8. Different colours correspond 
to different sites i and different averages fii. The Poissonian 
distributions (•) V^^ik) for a close system with density /ii, 
are compared to the curves (x), describing the open system 
with gradient; lines are guide to the eyes. The inset shows the 
average number (n) as a function of the site index i: (n) = 8 
approximately corresponds to i = 12. The largest discrepan- 
cies between P(i, k) and V^^ (k) occur for sites corresponding 
to occupation numbers around Nc- 

as pointed out in the inset). Otherwise stated, the dif- 
ference is appreciable in a region around ic- Moreover, 
experiments show that such a region sensibly shrinks as 
the size gets larger. 



C. Fluctuations 

Let us study the occupation number {ni) and its fluc- 
tuations F{ni) as functions of the system parameters 
{L,Ni,Nl,Nc). The interest in this point depends on 
the possibility, for interacting systems, of an increase in 
fluctuations induced by the presence of a gradient, as 
found e.g. in [ill, [l^. More precisely, the system consid- 
ered in [ll| was an Ising model on a cyhndrical lattice, 
in contact with two thermostats at temperatures Ti and 
T2. In the presence of a heat flow (Ti ^ T2), it was found 
that fluctuations (generically denoted as F) relevant to 
several observables all satisfy the following 

AF = Fflow — -Fno-flow > , 

where inflow and -Fko-flow represent the fluctuations in the 
system with and without flow, respectively. Moreover, 
such difference is especially important in a domain 5E of 
energies around Ec, the critical energy of the Ising system 
at equilibrium. 

Clearly, analogous comparisons can be carried out for 
our RW model, by considering fluctuations in the oc- 
cupation number for open steady systems in the absence 
(TVi = Nl) or presence (iVi ^ Nl) of a particle flow. No- 
tice that, apart from negligible border effects, the former 
case is well described by the close system of Section IIIII 



FIG. 4: Fluctuations in the local occupation number F(ni) 
as a function of (n^). Different sizes are depicted, a shown by 
the legend. iVi = 2, iVa = 10. 
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FIG. 5: Fluctuations in the local occupation number Fiui) as 
a function of (rii). Different gradients are depicted, as shown 
by the legend. L = 32. 

where the exact value of the ffuctuations can be evaluated 
from the Poissonian distribution (i.e. i^no-flow = i^)), al- 
lowing for an accurate insight into the problem. 

Figures [4] and [5] show AF as a function of {n) for dif- 
ferent choices of the system parameters. Actually, when 
a steady particle flow is established, («) AF > 0, {ii) at a 
density (n) ~ Nc deviations from the equilibrium values 
are stronger, and (iii) the overall effect is more signiflcant 
in a region 6N = {Nc — Si,Nc + 62)- The analogy with 
the Ising system is notable. 

Figures m and [H] also highlight the role of the size and of 
the density gradient: by enlarging L or by reducing Nl — 
Nl, the region SN gets smaller and the discrepancy AF 
shrinks. Again, this effect is consistent with the results 
found for the ferromagnetic model undergoing a steady 
heat flow. 

Thus, both the ferromagnetic and the RW systems ex- 
hibit an increase in fluctuations when a gradient is es- 
tablished. However, such a gradient is not a sufficient 
cause: sites have to interact so that a dynamical slowing 
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FIG. 6: AF for systems of different sizes and equal boundary 
conditions, i.e. A'^i = 2 and N2 — 10. The peak corresponds 
to Nc — 8. Different sizes are shown with different colours 
and symbols. Lines are guide for the eye. 



down takes place for a value of the parameter (energy and 
particle density respectively). For the Ising model, such 
a special value of the energy corresponds to the critical 
point. As for RW, it is the parameter Nc in the interac- 
tion defined by Eq. [TH 

Now, Ising and RW systems display an important dif- 
ference: while for L — > cx3 the former exhibits a true 
critical point, the latter is never critical, since Nc does 
not correspond to any singular behaviour. In particular, 
(to be related to the diffusivity) is always finite. 

We can conclude that the similar increase in fluctu- 
ations, observed in both systems, is not really due to 
criticality but only to a dynamical slowing down. This 
solves the apparent contrast remarked in the Ising case 
between the mesoscopic nature of the effect and the ther- 
modynamic nature of the critical point. 



D. Scaling 

Previous observations can be stressed by looking di- 
rectly at the scaling behaviour. In general, as well known, 
fluctuations are sensitive to criticality. The functional de- 
pendence on the size L, for instance, is different at the 
critical point. This happened indeed in the Ising model at 
Ec [m . Now, since in the present case Nc is not a genuine 
critical point, as L grows we expect, conversely, a regular 
behaviour for all occupation numbers, those approaching 
Nc included. Such an expectation is confirmed by exper- 
iments reported in Figure [HI where AF{{n)) goes to zero 
as L gets larger according to the scaling law: A_F ~ 1/L, 
with an excellent data collapse. The only relevant effect 
appearing at Nc regards the amplitude of A_F, which has 
a maximum, not a distinct scaling law. 



E. Spectral features 

Time series can provide further information on the 
pseudo criticality around Nc- It is worth summarizing 
how this kind of analysis is performed. 
For a given site i, consider the sequence of occupation 
numbers: 

{ni(ii),nj(i2), ...,nj(i„)}: 

where ti,t2, ■■■,tn are instants of time after the onset of 
a steady state. This definition may be implemented by 
introducing a time delay At, and the corresponding time 
series: 

{ni{ti),ni{ti + At),ni{ti + 2At), ...,n,{ti + kAt), ...}. 

Clearly, the delay or sampling parameter At plays an 
important role in experiments, since time correlations are 
strong in our model, and a sufficiently long interval is 
needed before a given configuration significantly changes. 

Then, by Fast Fourier Transform, from time series 
we get power spectra <S'(/), i.e. the square of absolute 
Fourier transform amplitudes in the frequency domain. 
The so-called colour exponent a describes the possible 
power-law decay S{f) ~ /" of the spectrum. It can be 
obtained as the angular coefficient of the linear fit in the 
log-log plot of S'(/) vs. /. 

Of course, a special value is a = associated to white 
noise, decorrelation, randomness. Values a < qualify 
coloured noises corresponding to different types of tem- 
poral correlations. In particular, a = —1 ("pink" or 1// 
noise) implies an extremely slow decay of correlations. 

As it is well known, a general dynamic theory of 
coloured noise is still lacking (see however [131 )■ The old 
suggestion of Van der Ziegle, getting 1// noise from the 
superposition of independent Poisson processes [2l|, , 
could be of some interest in our case. 

In simulations, for every At, exponents are evaluated 
and averaged over several runs starting from different 
initial conditions, up to stabilization. The dependence 
of the noise on interaction is shown in Fig. [71 where ex- 
ponents a versus (n) are plotted for analogous systems 
made up of interacting and non-interacting (i.e. isotrop- 
ically diffusing) particles; several choices of At are also 
considered. Systems of larger sizes, not reported here, 
display qualitatively alike outlines. 

The presence of interaction clearly shifts a towards 
smaller values, but only in a domain around Nc largely 
overlapping the previously considered domain 6N. In- 
deed, a is minimum for n ~ Nc- This confirms that 
around Nc correlations are stronger and the dynamics is 
slowed down. Moreover, the comparison with the non- 
interacting system allows to figure out the background 
time-correlation which, as expected, decreases as At is 
larger, up to a = 0. 

A less trivial point has to be underlined: the difference 
between interacting and non-interacting case is not a 
monotonic function of At and, therefore, neither of a 
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FIG. 7: Color exponent for a system of interacting (red) and 
non interacting (black) particles, both sized L = 16 and with 
Ni = 0.5, Nl = 20.0; for the former Nc = 8.0. Different time 
delays At are depicted, as shown in the legend. 

itself. Actually, there exists an optimal time delay At 
maximizing the effect of interaction over time correlation. 
Such a special delay At seems difficult to be interpreted, 
and could suggest that there exist further aspects, besides 
the superposition of Poisson processes, characterizing the 
noise in our model. 



VI. CONCLUSIONS AND PERSPECTIVES 

We have introduced a system of interacting random 
walkers modelling several situations where diffusivity de- 



pends on the density of the diffusing entities. Our model 
constitutes a special case of misanthrope processes, char- 
acterized by a totalistic hopping amplitude. We have 
shown analytically that the equilibrium distribution is 
factorized into a product of Poissonian functions. 

Despite a simple equilibrium behaviour, far from equi- 
librium the study of the system presents some non-trivial 
tasks. Indeed, in principle, local properties such as cor- 
relations and factorizability are deeply influenced by non 
equilibrium. In our case, we succesfuUy explored the reli- 
ability of a perturbative approach to the problem. Since, 
on one hand, the equilibrium properties are analytically 
known and, on the other, interaction still plays a rele- 
vant role, our model can be considered as an interesting 
benchmark for the study of non-equilibrium effects. By 
recovering analogous results in magnetic models, we have 
clarified some open problems, concluding, in particular, 
that a density dependent diffusivity and the presence of 
a dynamical slowing down are necessary premises to the 
amplification of fluctuations, while, on the contrary, the 
existence of a true critical point may cooperate but it is 
not strictly necessary. 

Finally, even if the passage to regular lattices of higher 
dimensionality does not present any substantial novelty, 
we evidenced a possible deep interplay between time scal- 
ing, substrate topology and local interactions. Conse- 
quently, there still exist in this context a number of 
geometro-dynamical aspects which have to be developed, 
understood and classified. The non trivial influence of in- 
teraction on the noise can be considered as a notable hint 
for the relevance of the matter. 
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